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TRANSLATION  AUTHOR'S  SUMMARY 

It  is  shown  that  there  is  an  interaction  between  pressure  gradients 
normal  and  parallel  to  the  surface  in  supersonic  flows  and  that  the  effect  of 
longitudinal  curvature  increases  strongly  with  Mach  number.  The  second  order 
boundary  layer  equations  and  the  matching  conditions  with  the  inviscid  flow 
are  discussed.  As  well  a*  longitudinal  curvature,  transverse  curvature  and 
variations  in  stagnation  pressure  and  enthalpy,  are  considered  as  second  order 
effects.  Integral  relationships  of  the  boundary  layer  equations  (momentum, 
energy  equation,  etc.)  are  given.  The  effects  are  demonstrated  by  a  numerical 
example  of  a  laminar  boundary  layer  profile  and  possibilities  for  an  approximate 
solution  of  the  equations  are  discussed. 
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1  IHTROPUCTION 

Methods  of  calculation  for  boundary  layers  have  attained  &  high  degree 
of  perfection.  The  main  questions  still  unsolved  are  those  in  which  simple 
boundary  layer  theory  is  not  sufficient. 

Prandtl's  boundary  layer  theory  Is  based  on  the  assumption  of  constant 
pressure  across  the  boundary  layer.  In  the  case  of  boundary  layers  on  curved 
surfaces  this  assumption  is  often  not  Justified.  In  recent  times  attempts 
have  been  made  to  determine  the  curvature  effect  using  second  order  boundary 
layer  theory  (see,  for  example,  Murphy1,  Yen  and  Toba2,  Hayasi^,  van  Dyke^'  ^ 
and  Massey  and  Clayton  ) .  Moat  of  the  work  is  limited  to  the  case  of 
incompressible  flow.  It  is  shown  that  significant  effects  may  be  expected 
when  the  thickness  of  the  boundary  layer  is  coop  arable  with  the  radius  of 
curvature  of  the  surface.  Much  stronger  effects  arise,  however,  in  flowa  with 
high  Mach  numbers. 

A  feature  of  compressible  flow  is  the  fact  that  pressure  gradients  normal 

to  the  surface  produced  by  centrifugal  forces  result  in  variations  of  the 

density  distribution.  These  directly  affect  the  force  equilibrium  end  thus 

influence  the  development  of  the  boundary  layer.  In  the  case  of  turbulent 

boundary  layers  the  turbulence  mechanism  is  also  influenced  in  a  complicated 

manner  so  that,  even  with  moderate  Mach  numbers,  strong  induced  effects  are 

observed,  as  has  been  shown  by  boundary  layer  studies  on  a  waistsd  body  of 

„  7 

revolution  which  were  carried  out  Jointly  by  R.A.E.  Bedford  and  AVA  Gottingen  . 

Further  flow  problems  exist  which  are  outside  the  scope  of  boundary  layer 
theory.  These  occur  primarily  in  flows  with  separation  and  near  the  trailing 
edge  of  bodies  (see  the  work  of  Ruchemann1 1 ) .  For  this  ■type  of  problem  a  theory 
of  order  higher  than  second  is  necessary  which  leads  to  differential  equations 
of  the  elliptic  type.  For  the  second-order  boundary  layer  theory  used  in  the 
present  investigation,  the  equations  form  a  problem  which  can  be  solved  by 
•marching*  downstream  as  in  the  case  of  the  simple  boundary  layer  theory.  The 
details  are  limited  to  laminar  boundary  layers,  although  many  are  also  valid 
for  turbulent  boundary  layers. 

2  THE  EFFECT  OF  SURFACE  CURVATURE 

In  order  to  emphasize  the  significance  of  the  surface  curvature.  Fig .1 
recalls  the  fact  that  for  supersonic  flow  around  a  body  there  is  a  causal 
relationship  between  the  pressure  gradients  normal  end  parallel  to  the  wall. 

The  pressure  variations  axe  propagated  along  Mach  lines.  According  to  simple 
theory  the  following  relationship  exists  between  the  pressure  gradients 


Tor  high  Mach  numbers  Ma  this  is  reduced  to 


!>-!!“»  •  (2) 


The  pressure  gradient  normal  to  the  surface  increases  proportionately  to 
the  pressure  gradient  in  the  flov  direction,  and  increases  strongly  with  Mach 
number  reaching  a  multiple  of  dp/ds. 

It  is  also  ohcum  that  the  relative  pressure  change  in  the  boundary  layer 
can  be  quite  considerable .  The  approximation  for  the  pressure  gradient  normal 
to  the  surface  gives 


.  „  2 
„  eji_ 
dy  R 


(5) 


In  the  case  of  an  ideal  gas,  the  integration  of  equation  (3)  gives,  for 
the  ratio  of  the  local  pressure  to  the  wall  pressure, 


where  7  is  the  ratio  of  specific  heats .  It  can  be  estimated  from  this  equation 
that,  for  exanple,  for  Ma  ■  10  the  pressure  at  the  boundary  of  a  flow  layer  of 
thickness  only  11&  of  Rw  is  approximately  4  times  the  surface  pressure. 

The  distribution  of  pressure  in  the  boundary  layer  on  a  convex  surface  is 
shown  in  Fig.2.  Since  the  actual  velocity  u  is  different  from  the  value  u* 
of  the  invlscid  flow,  different  pressures  are  obtained  for  the  two  flows.  For 
convex  surfaces  the  actual  pressure  is  greater  than  for  inviscid  flow.  The 
difference  can  be  considerable .  Consequently  a  boundary  layer  theory,  which 
neglects  tbs  influence  of  curvature  can  give,  under  some  circumstances, 
erroneous  results  at  high  Mach  numbers.  In  the  following  the  longitudinal 
curvature,  the  transverse  curvature,  the  variation  of  the  total  pressure 
(entropy)  end  total  enthalpy  are  considered  as  second-order  effects.  Slip  and 
tesperatuxe  Jung)  at  the  surface  are  not  considered,  although  inclusion  of  these 
effects  does  not  present  any  fundamental  difficulty. 
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3  BOUNDARY  LAYER  EQUATIONS 

The  general  formulation  of  the  second  order  boundary  layer  theory  corres¬ 
ponds  exactly  to  Prandtl’s  concept  (see  Schlichting®) :  friction  stresses  end 
heat  conduction  are  only  considered  in  a  relatively  thin  boundary  layer  near 
the  surface  of  the  body.  Inside  this  boundary  layer  certain  terms  of  the 
Navier-Stokes  differential  equations  and  the  energy  equation  may  be  neglected. 
In  the  remaining  flow  field  friction  stresses  and  heat  conduction  are 
neglected.  We  shall  discuss  the  boundary  layer  equations  for  steady  axial 
flow  past  a  body  of  revolution. 

An  orthogonal  axis  system  a,  y  is  introduced  on  a  curved,  axisymmetric 
surface  (Fig .3).  The  generator  of  the  surface  is  given  by  r^(x) .  The 
geometrical  relationships  are  specified  on  the  figure.  The  boundary  layer 
simplifications  derive  fretn  the  assumptions 

v «  u,  d/ds  «  d/dy  .  (5) 

Deviating  from  the  usual  boundary  layer  theory,  we  do  not  retain  the 
restrictions 


(6) 

where  6  is  the  boundary  layer  thickness.  The  boundary  layer  equations  for 
this  ease  have  already  been  treated  by  van  Dyke"* .  We  use  the  continuity, 
momentum  and  energy  equations  in  the  following  form: 

Continuity: 


V\  «  «/»■.  6  «  * 


affiurl  a  (p.vrR).  0 

ai  +  dy  u 


(7) 


Equations  of  motion: 
meridional: 


RRw 


d(p  u2  r)  d  (puvrR2) 
5^  •  '  dy 


d(tr  R2^ 
dy 


(6) 


normal  to  the  surface 


2 

p  u  dp 

R  "  -8? 


(9) 
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Energy: 


a(puhor)  3(pvhorR)  a(qrR)  afrurR) 
%  5s  +  dy  “  dy  dy 


(10) 


The  total  enthalpy  hQ  is  defined  by 


h 


o 


(11) 


The  shear  stress  in  the  laminar  case  is 


x 


(12) 


and  the  heat  flux  is  calculated  from 


2L  dh 

q  cpdy 


(13) 


where  p  is  the  viscosity,  X  the  thermal  conductivity,  the  specific 

heat  at  constant  pressure  and  h  the  enthalpy. 

Apart  from  the  occurrence  of  radii  of  curvature,  the  equations  differ 
from  those  of  simple  boundary  layer  theory  mainly  by  the  addition  of  the  equa¬ 
tion  for  the  change  of  pressure  normal  to  the  surface,  equation  (9) .  This 
raises  by  one  the  order  of  the  differential  equation  system.  The  mathematical 
consequences  for  the  behaviour  of  the  equation  system  have  not  yet  been 
coop  lately  investigated,  but  it  appears  that  the  fundamental  properties  remain 
unchanged  -  This  suggests  the  use  of  similar  methods  of  solution  as  in  usual 
boundary  layer  theory. 

Up  to  now  no  restrictions  have  been  made  on  the  behaviour  of  the  gas, 
apart  from  the  fact  that  it  is  taken  as  a  continuum.  At  best  the  expression 
for  q  according  to  equation  (13)  is  only  valid  in  the  case  of  thermodynamic 
and  chemical  equilibrium.  The  behaviour  of  the  gas  must  be  established  in 
order  to  solve  the  system  of  equations  and  for  the  sake  of  simplicity,  thermo¬ 
dynamic  and  chemical  equilibrium  will  be  assumed  here.  Pressure,  density  and 
enthalpy  are  then  related  by  an  equation  of  state  of  the  usual  form.  The 
a  inpleat  case  is  to  take  the  equation  of  state  of  an  ideal  gas;  however,  it 
ia  not  necessary  to  be  limited  to  this  assumption.  The  addition  of  further 
differential  equations  to  take  into  account  real  gas  effects  makes  the  equation 
system  quite  cccplicated  but  the  quoted  relationships  remain  valid. 
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4  BOUKDABY  C0HDITI0H3 

Uhder  the  assumption  of  an  impermeable  surface,  the  velocity  components 
must  fulfil  the  conditions  ■  0  and  vy  «  0  for  y  »  0.  Also  alternatively 
h  »  1^  or  q  «•  qw  can  be  given*. 

The  formulation  of  the  boundary  conditions  is  very  simple  with  first 

order  boundary  layer  theory.  Since  the  pressure  is  constant  through  the 

boundary  layer,  the  pressure  distribution  of  the  potential  flow  at 

the  wall  is  used  as  a  boundary  condition,  which  is  independent 

of  y.  in  the  case  of  the  second-order  theory  the  matching  to  the  inviacid 

flow  is  more  complicated  since  the  properties  axe  dependent  upon  y.  We  define 

a  thickness  6  such  that  friction  and  heat  conduction  terms  need  be  considered 
r 

only  for  0  <  y  <  and  are  negligible  for  y  *  6^  (Fig. 4).  Further,  we 
define  a  boundary  layer  thickness  8  >  6y  so  that  for  0  S  y  i  6  the  boundary 
layer  approximations  are  permissible.  For  y  *  d  the  complete  differential 
equations  of  the  in  viscid  flow  sire  used.  Thus,  in  the  overlapping  region 
6f  S  y  S  8,  boundary  layer  approximations  are  permissible  and  the  friction 
forces  end  heat  conduction  are  negligible.  The  boundary  conditions  for 
8r  <  y  S  6  then  reads 

u  ■  u*,  v  m  v*,  h  •  h*,  p  ■  p*  ,  (14) 

The  requirement  6y  i  8  defines  the  limit  for  the  applicability  of  the 
stated  theory,  i.e.  the  viscosity  must  be  low  enough  or  the  Reynolds  number 
high  enough  for  this  condition  to  be  fulfilled. 

The  in  viscid  flow  may  be  considered  as  continuing  up  to  the  surface  of 
the  body .  In  the  region  0  *  y  *  6  it  obeys  the  boundary  layer  equations  in 
which  the  friction  and  heat  conduction  terms  are  omitted.  It  is  advantageous 
to  describe  the  lnviscid  flow  by  the  tangential  velocity  u^J  or  the  pressure 
p*  or  h£  at  the  surface.  These  conditions,  however,  are  not  adequate  to 
specify  the  flow  field.  The  aquations  of  inviscid  flow  indicate  that  total 
pressure  p*  and  total  enthalpy  h*  are  constant  along  streamlines.  Thus 
p*  and  h*  must  be  known  u  functions  of  the  stream  function.  Only  in  special 
casas  (e.g.  irrotatlanal  flow)  are  these  parameters  constant. 


*  Taking  into  account  slip  and  tenperature  Jusp  at  the  surface,  uy  and  v 
are  given  from  a  linear  relationship  with  the  gradients  (du/dy)^,  (dh/8y)w 
and  (8h/8s)w  (see  Wuest^). 
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It  may  be  further  noted  that  the  requirement  v*  =  v  for  6f  S  y  S  t  ie 
not  compatible  with  the  condition  v*  =  0  for  y  =  0.  The  values  of  v*  can 
only  be  determined  in  the  case  of  a  known  boundary  layer.  In  many  instances 
the  displacement  effect  of  the  boundary  layer  on  the  external  flow  is  neglected. 
Possibilities  exist  of  introducing  the  displacement  effect  into  the  calculation, 
i.e.  in  an  iterative  correction  or  in  the  introduction  of  pressure  variations  in 
the  boundary  layer  equations  (interaction  theory);  this  will  not  be  discussed 
further  here . 

5  CTTEGBAL  EQUATIONS  OF  THE  BOUNDARY  LAYER  THEORY 

Among  the  mathematical  methods  for  the  solution  of  boundary  layer  equa¬ 
tions,  approximation  methods,  which  operate  with  integral  theorems  for  momentur^ 
energy  etc.,  are  of  particular  interest.  Such  methods  will  now  he  discussed. 

The  principle  of  the  integral  approach  is  that  we  consider  the  equilibrium 
of  the  momentum  and  energy  losses  produced  by  friction  forces  and  heat  conduc¬ 
tion  with  respect  to  the  inviacid  flow  considered  as  extended  to  the  wall.  In 
the  derivation  of  the  integral  conditions  we  assume  the  inviscid  flow  to  he 
known  and  require  that  in  the  region  0  i  y  i  6,  as  already  mentioned,  the 
boundary  layer  equations  are  adequate  without  the  friction  and  heat  conduction 
terms.  The  formal  derivation  is  mathematically  simple  and  need  not  he  explained 
further  here.  It  may  be  noted  that  in  the  boundary  layer  not  only  velocity, 
density  and  enthalpy  but  also  static  pressure  values  are  assumed  which  differ 
from  those  of  the  inviacid  flow.  After  rearranging  the  terms,  the  equations 
can  he  brought  into  forma  which  are  largely  similar  to  those  of  usual  boundary 
layer  theory,  and  which  reduce  to  these  if  the  radius  of  curvature  relative  to 
the  boundary  layer  thickness  tends  to  infinity. 

By  multiplication  of  the  equation  of  motion  by  um,  a  differential 
equation  for  a  loas  thickness  t> 2+m  can  be  established  for  each  value  of 
m.  For  m  =  0  this  gives  the  momentum  equation  in  the  usual  way;  for 
m  ■  1  we  have  the  integral  equation  for  mechanical  energy,  while  for  higher 
values  of  m  there  is  no  obvious  physical  explanation. 

Tbs  integral  form  of  the  boundary  layer  equation  is: 
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(20) 
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Ow  integral  parameter*  are  distinguished  from  thoee  of  ainjle  boundary 
layer  theory  (aee  Walt10:  only  by  th*  loaa  tl.-cknea.ee  6g+B  r  «nd  Sg+m>R- 
The  remaining  second  order  effects  are  alao  included  in  the  definition,  of  the 
loaa  thickness.  Cue  to  the  curvature  normal  to  the  streamlines,  the  radii 
ratio  r/r  appear  in  th*  integrand*  with  the  exception  of  ®g+m,r* 
case  of  6g+B  the  first  ten#  correapcoda  to  that  of  ainple  theory.  In  the 
second  teraTonly  properties  of  the  lnvi.cid  flow  occur.  This  term  disappears 
for  ieen tropic  flow  for  which  u*R  -  const.  In  the  third  term  the  difference 
between  the  static  pressures  appears.  The  expression  for  6^,,  in  the  inte¬ 
grands  of  which,  partial  differential  coefficients  also  arise,  is  written  with 
ft  in  a  common  term  like  that  of  usual  boundary  layer  theory.  For  m  -  0, 


*»• 


0. 


On  the  right  side  of  the  integral  equation  are  the  friction  terms .  The 
first  term  contains  th*  wall  shear  stress  and,  since  ^  -  0,  has  a  finite 
rain*  only  when  m  -  0  (momentum  equation);  in  this  case  the  second  term 
also  dis^pears.  Th*  latter  corresponds  to  the  dissipation  of  mechanical 
energy  for  m  •  1  • 

The  Integral  equation  of  total  energy  has  the  form: 


rsit'V:'."1 


dh  * 

(P  «>;  -g  6i 


(21) 


Where  q^  it  the  heat  flow  per  unit  are*  from  the  body*.  For  the  enthalpy 
loaa  thickness 


tha  first  term  in  the  square  brackets  again  correspond!  to  simple  boundary  layer 
theory.  The  second  tern  contains  only  properties  of  the  in  viscid  flow.  It  dis¬ 
appears  in  the  case  of  iaoenergatic  flow,  since  then  h*  ■  corst. 


*  If  we  relinquish  th*  no-slip  condition  at  the  surface  then  th*  term  +  xv\  i« 

added  to  the  right  hand  side  of  equation  (21)  and  the  first  term  of  the  right 
h.n*  tide  of  equation  (15)  has  a  finite  value  for  m  >  0.  The  equations  are 
then  valid  even  in  th*  case  of  a  temperature  Junp  at  the  surface . 
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Th#  lots  thicknesses  can  be  calculated  for  given  velocity  and  enthalpy 
profile# .  Aa  an  exaspla  we  will  show  calculated  reeulta  for  two-dimensional 
flow  over  convex,  plana  and  concave  (urfacea .  The  Pohlhausen  profile  assumed 
ia  eaaily  modified  to  natch  the  inviacid  flow 


The  equation#  of  atate  of  an  ideal  gaa  are  taken.  The  inviacid  flow  ia 
taken  aa  lrrotaticnal  with  hQ  «  conat .  The  velocity  profile#  deviate  only 
•lightly  from  each  otuer  for  a  relative  boundary  layer  thlckneea  of  6/R^  *  0.01 
(rig. 5).  On  the  other  hand  the  atatlc  preiaure  at  MgJ  •  8  undergo## 
considerable  variation.  The  marked  difference  between  actual  preaaure  and  the 
pressure  of  the  inviacid  flow  at  y  ■  0  is  particularly  noteworthy.  This 
difference  ia  shown  in  fig .6  aa  a  function  of  Mach  cumber. 

Fig. 7  shows  the  integrands  of  the  momentum  loss;  I21  is  the  momenta# 
loss  thickness  of  the  velocity,  1^  Is  the  contribution  of  the  statlo 
pressure.  This  fraction  has  a  large  value  in  the  example  shown. 

rig  .8  shows  the  thickness  ratio  4114  10041  coefficient  of 

friction  cf  (non dimen#  local  wall  shear  stress  »  2  0*«)  *)  00  4  logarithmic 

scala  as  a  function  of  the  Mach  number.  With  large  Mach  embers,  Ma£ 
appreciable  variations  due  to  curvature  of  the  well  are  shown. 

This  suggests  that  an  approximate  method  on  tbs  basis  of  a  ■  ingle - 
parameter  velocity  profile  should  he  sought,  where  the  Pohlhausen  polynmitl 
equation  can  ha  used.  A  combination  of  the  momentum  (m  ■  0)  and  mechanical 
energy  (a  >  1)  equations  appears  most  promising. 


Another  inportant  relationship  which  if  often  used  in  boundary  layer 
calculation  is  the  so-called  condition  of  compatibility  at  the  wall.  In  the 
case  of  second  order  boundary  layer  theory  this  read# 


...(24) 


Here  tbs  differential  coafflclent  of  an  integrated  boundary  layer 
parameter  appears,  so  that,  in  the  extension  of  the  Pohlhausen  method  to 
second  order  boundary  layer  theory,  e  first  order  differential  equation 
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trices  instead  of  an  algebraic  expression  for  the  shape  parameter.  Thus  in 
using  the  wall  compatibility  condition  and  dispensing  vitu  the  equation  for 
mechanical  energy;  a  system  of  two  simultaneous  differential  equations  is 
given. 

With  reference  to  tbs  temperature  profile,  the  aasuaption  that  the  same 
distribution  of  total  enthalpy  occurs  in  the  boundary  layer  as  in  the  inviscid 
flow  may  be  an  acceptable  approximation  in  the  case  of  insulated  surfaces 
(l.e.  h  ■  constant  for  isoenergetic  inviscid  flow).  In  the  case  of  heat 
transfer  the  integral  condition  for  the  energy  conservation,  equation  (Si), 
must  be  taken  as.  third  differential  equation  and  a  suitable  expression  for  the 
enthalpy  profile  Introduced. 

6  OOlKLOBIOff 

The  ispartance  of  second  order  boundary  layer  theory  in  hypersonic  flow 
ti..  been  indicated  and  an  example  given.  The  longitudinal  and  transverse 
curvature  of  the  surface  end  the  variations  of  tbs  static  pressure  and  static 
enthalpy  of  the  inviscid  solution  have  been  considered  as  second  order  effects . 
The  affect  of  the  longitudinal  curvature  has  been  specially  considered.  The 
mathematical  problem  has  bean  formulated  and  possibilities  of  approximate 
calculation  discussed.  In  order  to  have  a  basis  of  coup  arisen  for  the  calculated 
results,  more  accurate  methods  axe  necessary  based  either  on  a  larger  number  of 
integral  conditions  and  therefore  making  possible  a  larger  variation  of  the 
boundary  layer  profile  or  on  difference  methods  for  the  partial  differential 
equations.  The  quoted  integral  conditions  are  also  valid  for  turbulent 
boundary  layers;  however,  farther  data  on  the  velocity  profile,  friction 
coefficient  and  shear  stress  distribution  end  also  beat  transfer  coefficients 
ere  necessary  for  their  application. 
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